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ABSTRACT 

Rich clusters of galaxies are the most massive virialized systems known. Even 
though they contain only a small fraction of all galaxies, rich clusters provide a 
powerful tool for the study of galaxy formation, dark matter, large-scale struc- 
ture, and cosmology. 

Superclusters, the largest known systems of galaxies, extend to ~ 100/i _1 
Mpc in size and highlight the large-scale structure of the universe. This large- 
scale structure reflects initial conditions in the early universe and places strong 
contraints on models of galaxy formation and on cosmology. 

Some of the questions that can be addressed with clusters and superclusters of 
galaxies include: How did galaxies and larger structures form and evolve? What 
is the amount, composition, and distribution of matter in clusters and larger 
structures? How does the cluster mass density relate to the matter density in 
the universe? What constraints can the cluster and supercluster data place on 
cosmology? 

I will discuss in these lectures some of the properties of clusters and super- 
clusters of galaxies that can be used to investigate these topics. 

1. Introduction 

Clusters and superclusters of galaxies have been studied extensively both for their in- 
trinsic properties and to investigate the dark matter in the universe, the baryon content 
of the universe, large-scale structure, evolution, and cosmology. For previous reviews see 
Zwicky (1958), Bahcall (1977, 1988, 1996), Oort (1983), Dressier (1984), Rood (1988), and 
Peebles (1993). 

In these lectures I discuss the following topics and their implications for structure for- 
mation and cosmology. 

Section 2: Optical properties of galaxy clusters 
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A Hubble constant of H Q = 100/t km s 1 Mpc 1 is used throughout. 

2. Optical Properties of Galaxy Clusters 

2.1. Typical Properties of Clusters and Groups 

Clusters of galaxies are bound, virialized, high overdensity systems of galaxies, held 
together by the clusters self gravity. Rich clusters contain, by traditional definition (Abell 
1958), at least 30 galaxies brighter than m 3 + 2 m (where is the magnitude of the third 
brightest cluster member) within a radius of R ~ 1.5/i _1 Mpc of the cluster center. This 
galaxy count is generally defined as the richness of the cluster. The galaxies in rich clusters 
move with random peculiar velocities of typically ~ 750 km s _1 (median line-of-sight velocity 
dispersion). This motion corresponds to a typical rich cluster mass (within 1.5/i _1 Mpc) 
of ~ 5 x 10 14 /i _1 Mq. In addition to galaxies, all rich clusters contain an intracluster 
medium of hot plasma, extending as far as the main galaxy concentration (R ~ 1 .5/i —1 Mpc). 
The typical temperature of the hot intracluster gas is ~ 5 kev, with a range from ~ 2 to 
14 kev; the central gas density is ~ 10~ 3 electrons cm -3 . The hot plasma is detected 
through the luminous X-ray emission it produces by thermal bremsstrahlung radiation, with 
L x ~ 10 44 erg s _1 . 
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Like mountain peaks on earth, the high density rich clusters are relatively "rare" objects; 
they exhibit a spatial number density of ~ 10~ 5 clusters Mpc -3 , as compared with ~ 10~ 2 
galaxies Mpc -3 for the density of bright galaxies. 

The main properties of clusters and groups of galaxies are summarized in Table |l] (Bah- 
call 1996). The table lists the typical range and/or median value of each observed property. 
Groups and poor clusters, whose properties are also listed, provide a natural and continuous 
extension to lower richness, mass, size, and luminosity from the rich and rare clusters. 

In the following subsections I discuss in more detail some of these intrinsic cluster 
properties. 

2.2. Distribution of Clusters with Richness and Distance 

An illustration of the distribution of rich clusters with richness and distance is sum- 
marized in Table ^|. The table refers to the statistical sample of the Abell (1958) cluster 
catalog (i.e., richness class R > 1, corresponding to a richness threshold count of Nr > 50 
galaxies within 1.5/t -1 Mpc radius and magnitude m < m 3 + 2 m (see Table |l]); redshift range 
z ~ 0.02 to 0.2; and sky coverage 5 > —27° and |6| ;> 30°). The Abell catalog covers ~ 1/3 
of the entire sky to z ^ 0.2. Recent smaller automated surveys from digitized plates are 
reported by the Edinburgh-Durham catalog (EDCC; Lumsden et al. 1992), and the APM 
survey (Dalton et al. 1992). A large automated cluster catalog will be available in the near 
future from the Sloan Digital Sky Survey (SDSS). This and other planned surveys will allow 
a more accurate determination of the distribution of rich clusters with richness and distance 
and other statistical studies of clusters. 



2.3. Number Density of Clusters 

The number density of clusters is a strong function of cluster richness. Integrated cluster 
densities, n c (> Nr), representing the number density of clusters above a given richness 
threshold, and the associated mean cluster separation, d (= n~ 1//3 ), are listed in Table |] 
(Bahcall and Cen 1993). 
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Table 1: Typical Properties of Clusters and Groups. 



Property 


Rich clusters 


Croups 
and poor clusters 


Richness 6 


30-300 galaxies 


3-30 galaxies 


Radius c 


(1-2) h~ l Mpc 


(0.1-1) h' 1 Mpc 


Radial velocity dispersion 01 


400-1400 kins" 1 


100-500 kms" 1 


Radial velocity dispersion^ 


~ 750 kms" 1 


~ 250 kms" 1 


(median) 






Mass (r < 1.5/r 1 Mpc) e 


(lO^xlO 15 )/?- 1 M Q 


(lO^-lO 14 )/! -1 M 


Luminosity (B)f 


(6 x 10 n -6xl0 12 )/i~ 2 L G 


(10 10 - 5 -10 12 )/i- 2 L Q 


(r < 1.5k- 1 Mpc) 








_ , Qnn/i ka It 

~ oUU/2 JVIqI Liq 


„ , onnA A/< it 
~ zUU/i JVIqI Liq 


X-ray temperature' 4 


2-14 keV 


£ 2 keV 


X-ray luminosity^ 1 


(10 42 - 5 -10 45 )/i- 2 ergs" 1 


£ 10 43 /i" 2 ergs" 1 


Cluster number density* 


(10" 5 -10- 6 )/i 3 Mpc~ 3 


(10- 3 -10~ 5 )/i 3 Mpc 


Cluster correlation scale- 7 


{22±A)h- 1 Mpc (i2 > 1) 


(13±2)/i- 1 Mpc 


Fraction of galaxies in 


-5% 


~ 55% 


clusters or groups fc 







a In most entries, the typical range in the listed property or the 
median value is given. Groups and poor clusters are a natural and 
continuous extension to lower richness, mass, size, and luminosity 
from the rich and rare clusters. 

b Cluster richness: the number of cluster galaxies brighter than m 3 
+ 2 m (where is the magnitude of the third brightest cluster 
galaxy), and located within a 1.5/jT 1 Mpc radius of the cluster 
center (§gg). 

c The radius of the main concentration of galaxies (where, typically, 
the galaxy surface density drops to ~ 1% of the central core 
density). Many clusters and groups are embedded in larger scale 
structures (to tens of Mpc). 

d Typical observed range and median value for the radial (line-of- 
sight) velocity dispersion in groups and clusters (§ [2.9|) . 

e Dynamical mass range of clusters within 1.5/z" 1 Mpc radius 

rarmn. 



- 5 - 



Table 2: Distribution of Abell Clusters with Distance and 


Richness." 








Distance distribution 




Richness distribution 


D ( Zest ) N cl (R 


>1) 


R 


Nr N cl 


1 0.0283 


9 


(Of 


(30-49) (~ 10 3 ) 


2 0.0400 


2 


1 


50-79 1224 


3 0.0577 


33 


2 


80-129 383 


4 0.0787 


60 


3 


130-199 68 


5 0.131 


657 


4 


200-299 6 


6 0.198 


921 


5 


> 300 1 


Total 


1682 




Total (R>1) 1682 


Nearby redshift sample c ' d 


Distant projected sample^ 


D < 4 






D = 5 + 6 


A^i (total) 


104 




1574 


N cl (b > 30°) 


71 




984 


N cl (b < -30°) 


33 




563 


N cl (R = 1) 


82 




1125 


N cl (R > 2) 


22 




422 



a Statistical sample. |6| boundaries as given in Table 
1 of Abell (1958). Notation: D = distance group 
(defined by the estimated redshifts of the clusters); 
(z es t) = average estimated redshift from the mag- 
nitude of the tenth brightest galaxy; N d = number 
of clusters; R = richness class; N R = number of 
galaxies brighter than m 3 + 2 m within Ra = 1.5/i -1 
Mpc (richness count). 

b R = clusters are not part of the statistical sample. 

c Redshifts by Hoessel et al. (1980). 

d This sample is limited to |6| > 30° in addition to 
the |6| boundaries of the statistical sample. 
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2.4. Fraction of Galaxies in Clusters 

The fraction of galaxies in rich R > clusters is ~ 5% (within the Abell radius Ra = 
1 Mpc). The fraction of all galaxies that belong in clusters increases with increasing 
radius Ra and with decreasing cluster richness threshold. 

The average number of galaxies per cluster for R > clusters within 1.5/i -1 Mpc radius 
and m < m 3 + 2 m is (N R ) mcdimi ~ 50, or (N R ) mcan ~ 56. For R > 1 clusters the average is 
{Nr) median — 60, or {N R ) mcmi ~ 75. The number of galaxies increases to fainter luminosities 
following the Schechter (1976) luminosity function. 



2.5. Galaxy Overdensity in Rich Clusters 

The average number density of bright (> L*) galaxies in R > clusters (within Ra = 
1.5k- 1 Mpc) is 

n 9 (cluster) ~ 3h 3 galaxies Mpc -3 . (1) 
The average overall (field) number density of bright (^ L*) galaxies is 

n s (field) ~ 1.5 x 10 -2 /i 3 galaxies Mpc -3 . (2) 

The average galaxy overdensity in rich (R > 0) clusters (within 1.5/i _1 Mpc radius) is thus 

n g (cluster) /n g (field) ~ 200. (3) 

The typical threshold galaxy overdensity in clusters (within 1.5/i -1 Mpc radius) is 

R > clusters : n g (cluster) jn g (field) ^ 100, (4) 



Table 3: Number Density of Clusters. 



R 


Nr 


n c (> N R )h 3 (Mpc- 3 ) a 


d (> N R )h- 1 (Mpc) 


> o 


> 30 


13.5 xlO -6 


42 


> 1 


> 50 


6.0 xlO" 6 


55 


> 2 


> 80 


1.2 xlO -6 


94 


> 3 


> 130 


1.5 xl0~ 7 


188 



a Approximate uncertainties are 10 ±0 2 for the R > 0,1,2 
densities and 10 ±0 3 for R > 3. 
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R>1 clusters : n g (cluster) jn g (field) ^ 200. (5) 

The galaxy overdensity increases at smaller radii from the cluster center. The galaxy over- 
density in the cores of typical compact rich clusters is approximately 

nj(cluster core) fn g (field) ~ 10 4 -10 5 . (6) 



2.6. Density Profile 

The radial density distribution of galaxies in a rich cluster can be approximated by a 
bounded Emden isothermal profile (Zwicky 1957; Bahcall 1977), or by its King approximation 
(King 1972) in the central regions. 

In the central regions, the King approximation for the galaxy distribution is 

n g (r) = n° g (l + r 2 / R 2 C )~ 3/2 , spatial profile (7) 

S g (r) = S° g (l + r 2 /R 2 c )-\ projected profile. (8) 

n g (r) and S g (r) are, respectively, the space and projected profiles (of the number density of 
galaxies, or brightness), n g and S g are the respective central densities, and R c is the cluster 
core radius [where S(R C ) = S°/2]. Typical central densities and core radii of clusters are 
listed in the following subsection. The projected and space central densities relate as 

S° g = 2R c n° g . (9) 

A bounded Emden isothermal profile of galaxies in clusters yields a profile slope that varies 
approximately as (Bahcall 1977) 

S g (r ^ .R c /3) ~ constant, (10) 
S 9 (R C < r < 10R C ) oc r~ L6 ; (11) 

therefore 

n g (R c & r & 10R C ) oc r" 2 - 6 . (12) 

The galaxy-cluster cross-correlation function (Peebles 1980; Lilje and Efstathiou 1988) also 
represents the average radial density distribution of galaxies around clusters. For R > 1 
clusters, and r in h~ l Mpc, these references suggest, respectively 

£ 9C (r) ~ 130r~ 2 - 5 + 70r~ L7 (13) 

or 

-2.2 



e 9C (r) ~ 120r" 2 - 2 . (14) 
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The average density distribution profile of galaxies in clusters thus follows, approximately, 

n g {r) oc r~ 2A±0 - 2 (spatial), r > R c (15) 

S g (r) oc r - 1A±0 - 2 (projected), r > R c . (16) 

Some substructure (subclumping) in the distribution of galaxies exists in a significant fraction 
of rich clusters (~ 40%) (Geller 1990). 

2.7. Central Density and Core Size 

Central number density of galaxies in rich compact clusters (Bahcall 1975, 1977; Dressier 
1978) is (for galaxies in the brightest 3 magnitude range) 

n° g (Am ~ 3 m ) ~ 10 3 /i 3 galaxies Mpc -3 . (17) 

The central density reaches ~ 10 4 /i 3 galaxies Mpc -3 for the richest compact clusters. The 
typical central mass density in rich compact clusters, determined from cluster dynamics is 

Po(mass) ~ 9a 2 j4nGR 2 c (18) 
~ 4 x 1O 15 .M Mpc _3 [(<7 riC /10 3 kms~ 1 )/(i? c /0.2 Mpc)] 2 /i 2 

where a rtC is the radial (line-of-sight) central cluster velocity dispersion (in kms -1 ) and R c 
is the cluster core radius (in Mpc). 

Core radii of typical rich compact clusters, determined from the galaxy distribution 
(Bahcall 1975; Dressier 1978; Sarazin 1986) are in the range 

R c ~ (0.1 -0.25)/^ Mpc (19) 

Core radii of the X-ray emitting intracluster gas (|@ are 

,R C (X - rays) ~ (0.1 - O.S)/^ 1 Mpc (20) 

The core radius of the mass distribution determined from gravitational lensing observations 
of some clusters may be smaller, R c <, 50 kpc, than determined by the galaxy and gas 
distribution. 

The typical central density of the hot intracluster gas in rich clusters (§§) is 

n e ~ 10 -3 electrons cm -3 . (21) 
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2.8. Galactic Content in Rich Clusters 

The fraction of elliptical, SO, and spiral galaxies in rich clusters differs from that in the 
field, and depends on the classification type, or density, of the cluster (see § |2.11| ) (Bahcall 
1977; Oemler 1974; Dressier 1980). See Table |. 

The fraction of elliptical (E) and SO galaxies increases and the fraction of spirals de- 
creases toward the central cores of rich compact clusters. The fraction of spiral galaxies in 
the dense cores of some rich clusters (e.g., the Coma cluster) may be close to zero. 

The galactic content of clusters as represented in Table ^ is part of the general density- 
morphology relation of galaxies (Dressier 1980; Postman and Geller 1984); as the local 
density of galaxies increases, the fraction of E and SO galaxies increases and the fraction of 
spirals decreases. For local galaxy densities n g <, 5 galaxies Mpc -3 , the fractions remain 
approximately constant at the average "Field" fractions listed above. 



2.9. Velocity Dispersion 

The typical radial (line-of-sight) velocity dispersion of galaxies in rich clusters (median 
value) is 

<r r ~ 750 kms -1 . (22) 
The typical range of radial velocity dispersion in rich clusters (Struble and Rood 1991) is 

a r ~ 400-1400 km s _1 . (23) 

A weak correlation between a r and richness exists; richer clusters exhibit, on average, larger 
velocity dispersion (Bahcall 1981). The observed velocity dispersion of galaxies in rich clus- 
ters is generally consistent with the velocity implied by the X-ray temperature of the hot 



Table 4: Typical Galactic Content of Clusters (r < 1.5h 1 Mpc). 



Cluster type 


E 


SO 


Sp 


(E+S0)/Sp 


Regular clusters (cD) 


35% 


45% 


20% 


4.0 


Intermediate clusters (spiral-poor) 


20% 


50% 


30% 


2.3 


Irregular clusters (spiral-rich) 


15% 


35% 


50% 


1.0 


Field 


10% 


20% 


70% 


0.5 
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intracluster gas ( §3.5| ), as well as with the cluster velocity dispersion implied from obser- 
vations of gravitational lensing by clusters (except possibly in the central core). Velocity 
dispersion and temperature profiles as a function of distance from the cluster center have 
been measured only for a small number of clusters so far. The profiles are typically isothermal 
[of (r) ~ T x (r) ~ constant] for r < 0.5 — lh^ 1 Mpc, and drop somewhat at larger distances. 



2.10. Mass, Luminosity, and Mass-to-Luminosity Ratio 

The typical dynamical mass of rich clusters within 1.5/i _1 Mpc radius sphere (determined 
from the virial theorem for an isothermal distribution) is 

2^d.5ft-' Mpc) f^_\ 2 

v ~ ' G V 1000 / 

~ 0.4 x lO 15 ^" 1 ^® (for o r ~ 750 kms -1 ). 

The approximate range of masses for R > clusters (within 1.5/i -1 Mpc) is 

M d (< 1.5) ~ (0.1-2) x lO 15 /^ 1 M Q . (25) 

Comparable masses are obtained using the X-ray temperature and distribution of the hot 
intracluster gas (Hughes 1989; Bahcall and Cen 1993; Lubin and Bahcall 1993; §|3|). 

The typical (median) blue luminosity of rich clusters (within 1.5/t -1 Mpc) is 

L d (< 1.5) ~ 10 12 /i~ 2 L Q (26) 
The approximate range of rich cluster blue luminosities is 

L c i(< 1.5) ~ (0.6-6) x 10 12 h- 2 L Q (27) 
The typical mass-to-luminosity ratio of rich clusters is thus 

{M/L B ) cl ~ 300/i (M Q /L Q ). (28) 
The inferred mass-density in the universe based on cluster dynamics is 

fidyn ~ 0.2 (29) 

(if mass follows light, M oc L, on scales lh^ 1 Mpc). Q = 1 corresponds to the critical 
mass-density needed for a closed universe and M./L B (ft = 1) ~ 1500/i. 
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2.11. Cluster Classification 

Rich clusters are classified in a sequence ranging from early- to late-type clusters, or 
equivalently, from regular to irregular clusters. Many cluster properties (shape, concentra- 
tion, dominance of brightest galaxy, galactic content, density profile, and radio and X-ray 
emission) are correlated with position in this sequence. A summary of the sequence and 
its related properties is given in Table |^. Some specific classification systems include the 
Bautz-Morgan (BM) System (Bautz and Morgan 1970), which classifies clusters based on 
the relative contrast (dominance in extent and brightness) of the brightest galaxy to the 
other galaxies in the cluster, ranging from type I to III in decreasing order of dominance; 
and the Rood-Sastry (RS) system (Rood and Sastry 1971) which classifies clusters based on 
the distribution of the ten brightest members (from cD, to binary (B), core (C), line(L), flat 
(F), and irregular (I)). 

3. X-Ray Properties of Galaxy Clusters 

3.1. X-Ray Emission from Clusters 

All rich clusters of galaxies produce extended X-ray emission due to thermal bremsstrahlung 
radiation from a hot intracluster gas (Jones and Forman 1984; Sarazin 1986; David et al. 
1993; Edge et al. 1990, 1991; Henry et al. 1991, 1992; Burg et al. 1994). The cluster X-ray 
luminosity emitted in the photon energy band E\ to E2 by thermal bremsstrahlung from a 
hot (T x degrees) intracluster gas of uniform electron density n e and a radius i? x is 

L x oc n 2 e R 3 X- 5 g(e~ E ^ - e~ E ^). (30) 

The Gaunt factor correction g (of order unity) is a slowly varying function of temperature 
and energy. The bolometric thermal bremsstrahlung luminosity of a cluster core can be 
approximated by 

L x (core) ~ 1.4 x 10 42 n e (cm- 3 ) 2 J R c (kpc) 3 A;T x (keV) a5 /i- 2 ergs" 1 . (31) 
Some of the main properties of the hot intracluster gas are summarized below. 

3.2. X-Ray Properties of Clusters 

Some of the main properties of the X-ray emission from rich clusters of galaxies are 
summarized in Table ^ 
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Table 5: Cluster Classification and Related Characteristics. 





Regular (Early) 


Intermediate 


Irregular (late) 


Property 


type clusters 


clusters 


type clusters 


Zwicky type 


Compact 


Medium-compact 


Open 


BM type 


I, I— II, II 


(II), II III 


(II III), III 


RS type 


cD,B,(L,C) 


(L),(F),(C) 


(F),I 


Shape symmetry 


Symmetrical 


Intermediate 


Irregular shape 


Central concentration 


High 


Moderate 


Low 


Galactic content 


Elliptical-rich 


Spiral-poor 


Spiral-rich 


E fraction 


35% 


20% 


15% 


SO fraction 


45% 


50% 


35% 


Sp fraction 


20% 


30% 


50% 


E:S0:Sp 


3:4:2 


2:5:3 


1:2:3 


Radio emission 


~ 50% detection 


~ 50% detection 


~ 25% detection 


X-ray luminosity 


High 


Intermediate 


Low 


Fraction of clusters 


~ 1/3 


~ 1/3 


~ 1/3 


Examples 


A401, Coma 


A194 


Virgo, A1228 
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Table 6: X-Ray Properties of Rich Clusters. 



Property 


Typical value or range 


Notes 


L x (2-10 keV) 


~ (10 42 - 5 -10 45 )/i~ 2 ergs" 1 


a 




J x (r) oc [1 + (r/R c ) 2 ]- W+1/2 


b 


(P) 


~ 0.7 


c 


Pgas(r) 


p g as(r) oc [1 + (r/R c )*]- 3 ^ 






oc [1 + (r/i?,) 2 ]- 1 


d 


kT x 


~ 2-14 keV 


e 


T x 


~ 2 x 10 7 -10 8 K 


e 


Aspect — fcT x /Atm p 


~ 1 


f 


R c (x) 


~ (O.l-O.S)/^ 1 Mpc 


9 


n e 


~ 3 x lO^h 1 / 2 cm' 3 


h 


M gas (£ 1.5/i- 1 Mpc) 


~ 1O 13 - 5 A^ [range: (1O 13 -1O 14 )A1 


h -2.5] i 


M gas /M cl (£ 1.5/r 1 Mpc) 


~ 0.07 (range: 0.03-0.15 h" 1 - 5 ) 


i 


Iron abundance 


~ 0.3 solar (range: 0.2-0.5) 


3 



a The X-ray luminosity of clusters (2-10 keV band). (L x ) increases 
with cluster richness and with cluster type (toward compact, 
elliptical-rich clusters) (Bahcall 1977a,b; Sarazin 1986; Edge et 
al. 1991; Jones and Forman 1992; David et al. 1993; Burg et al. 
1994). 

b X-ray surface brightness distribution, J x (r); R c is the cluster core 
radius. 

c Mean ((3) from observations of X-ray brightness profiles (Jones 
and Forman 1984; Sarazin 1986). 

d Implied spatial density profile of the hot gas in the cluster [from 
I x (r); isothermal]. 

e Range of observed X-ray gas temperature in rich clusters (Edge 
et al. 1990; Henry and Arnaud 1991; Arnaud et al. 1992). 

f Aspect is the ratio of galaxy to gas velocity dispersion: /j, is mean 
molecular weight in amu (fi ~ 0.6), m p is mass of the proton, av- 
is radial velocity dispersion of galaxies in the cluster, and T x is 
the X-ray temperature of the gas (Lubin and Bahcall 1993). 

9 Cluster core radius determined from the X-ray distribution in the 
cluster (Jones and Forman 1992). 

h Typical intracluster gas density in rich cluster cores (Jones and 
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3.3. The Intracluster Gas: Some Relevant Questions 

Some of the fundamental questions about the intracluster gas relate to its origin, evo- 
lution, metal enrichment, hydrodynamical state in the cluster, and its relation to the distri- 
bution of galaxies and mass in the cluster. I list below some of the relevant questions and 
discuss a selection of these topics in the following subsections. Some of the questions posed 
do not yet have sufficient observational constraints to suggest a possible solution. 

• What is the hydrodynamical state of the hot gas in clusters? Is it in approximate 
hydrostatic equilibrium with the cluster potential? 

• What is the relation of the intracluster gas to the galaxies and mass in the cluster? 
For example (the subscripts below refer to gas, galaxies, and mass, respectively): 

— Density profiles: p gas (r) vs. p ga \(r) vs. p m (r)7 

— Temperature-velocity relation: T vs. a ga i vs. cr m ? 

— Mass: M gas vs. M gal vs. M d ? 

— Profiles of the above properties: 

T(r) vs. Og a i(r) vs. cr m (r)? 

M gas (r) vs. M gal (r) vs. M cl (r)? 

— Core radii: _R c (gas) vs. i? c (gal) vs. R c (m)l 

— Luminosity versus galaxy-type relation: L x vs. spiral fraction? 

• What is the origin of the hot intracluster gas? 

• What is the origin of the metal enrichment of the gas? 

• What is the evolution of the intracluster gas? 

3.4. The Intracluster Gas: Hydrostatic Equilibrium? 

The standard model of clusters assumes that both the gas and the galaxies are in 
approximate hydrostatic equilibrium with the binding cluster potential (Bahcall and Sarazin 
1977; Forman and Jones 1984; Sarazin 1986; Evrard 1990; Bahcall and Lubin 1994). In this 
model the gas distribution obeys 

dP gas = GM cl (< r)p, 
dr r 2 



(32) 
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where P gas and p gas are the gas pressure and density, and M d (< r) is the total cluster binding 
mass within a radius r. The cluster mass can thus be represented as 

Md (< r) = f + r , (33) 

pm p G V amr amr y 

where T is the gas temperature and /zm p is the mean particle mass of the gas. 

The galaxies in the cluster respond to the same gravitational field, and they satisfy 

-^>= 4(^+^+4* (34) 

where cr T is the radial velocity dispersion of galaxies in the cluster, p m \{r) is the galaxy 
density profile, and A represents a possible anisotropy in the galaxy velocity distribution 
[A = 1 — (a t /a r ) 2 , where t and r represent the tangential and radial velocity components]. 

The above two relations yield 

p _ cr 2 dlnp gas (r) I dlnr + dlnT / dlnr 

spec kT / nvrip dlnp ga \(r) / dlnr + dlna 2 / dlnr + 2A 

where the /3 spec parameter, defined by the left side of the above relation, can be determined 
directly from observations of cluster velocity dispersions and gas temperatures. The P spec 
parameter represents the ratio of energy per unit mass in the galaxies to that in the gas. 
Observations of a large sample of clusters yield a mean best-fit value of /3 spec ~ 1 ± 0.1 
(Lubin and Bahcall 1993; see also §|3.5|). This suggests that, on average, the gas and galaxies 
follow each other with comparable energies (of ~ kT/pm p ). The observed mean value 
Aspec — 1 =t 0.1 is consistent with the value of j3 determined from the right-hand side of 
the P relation (referred to as /3a t , and determined from the gas and galaxy density profile 
fits). Using p gas (r) oc r~ 2 (jQ and p gal (r) oc r ~ 2A±0 - 2 (Q), one finds p m ~ 0.85 ± 0.1 
(for an isothermal distribution) (Bahcall and Lubin 1994). The above consistency supports 
the assumption that the gas is in an approximate hydrostatic equilibrium with the cluster 
potential, and suggests that the galaxies and gas approximately trace each other in the 
clusters. 



3.5. The Relation between Gas and Galaxies 



The hot intracluster gas in rich clusters appears to trace reasonably well the galaxies in 
the clusters, and — with larger uncertainty — also the cluster mass. 
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kT (keV) 



Fig. 1. — Cluster radial velocity dispersion (ay) vs. gas temperature (kT) for 41 clusters 
(Lubin and Bahcall 1993). The best-fit (3 = a^/(kT/fxm p ) lines are shown by the solid and 
dotted curves, with f3 ~ 1. The /3 ~ 0.5 line previously proposed for a velocity bias in 
clusters is shown by the dashed curve; the velocity bias is inconsistent with the data. 
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Velocity-Temperature relation The galaxy velocity dispersion in clusters is well corre- 
lated with the temperature of the intracluster gas; it is observed (Fig. [[]) that a 2 ~ kT j \im p 
(Lubin and Bahcall 1993). The best-fit a—T relation is listed in § |3.7| . The observed correla- 
tion indicates that, on average, the energy per unit mass in the gas and in the galaxies is the 
same. Figure [l] shows that, unlike previous expectations, the galaxy velocities (and therefore 
the implied cluster mass) are not biased low with respect to the gas (and, by indirect impli- 



cations, with respect to the cluster mass; see also § |3.4|) . Results from gravitational lensing 
by clusters also suggest that no significant velocity bias exists in clusters, and that the gas, 
galaxies, and mass provide consistent tracers of the clusters. Cosmological simulations of 
clusters (Lubin et al. 1996) produce a-T correlations that match well the data in Figure [L| 



Density Profiles The gas density profile in clusters follows 

Pgas(r)^p gas (o)[l + (r/ J R c ) 2 l" 1 , (36) 



with core radii in the range R c ~ 0.1-0.3/i 1 Mpc (§3.2). This implies p gas ( r ) oc r 2 for 
R c <r < 1.5/T 1 Mpc. 



The galaxy density profile in clusters follows approximately 

Pg a i(r) oc r" 2 - 4±a2 R c <r^ 1.5/T 1 Mpc (37) 

with core radii R c ~ 0.1 - 0.25/T 1 Mpc (§[T7J). 

The mass density profile in clusters is less well established, but initial results from 
gravitational lensing distortions of background galaxies by foreground clusters suggest that 
the mass profile is consistent with the galaxy density profile (Tyson and Frische 1996). In 
the small central core regions of some clusters (r <^ 100 kpc), the mass distribution may be 
more compact than the gas or galaxies, with a small mass core radius of R c (m) < 50/T 1 kpc. 
The results for the overall cluster, however, suggest that the distributions of gas, galaxies, 
and mass are similar (with the gas distribution possibly somewhat more extended than the 
galaxies, as seen by the mean density slopes above). 



Beta-Discrepancy The mean /3 spec = a 2 / 'kT '/ [im v ~ 1 result discussed above, combined 
with the similarity of the gas and galaxy density profile slopes (that yields /3st ~ 0.85 ± 0.1; 
§ |3~5| ) show that the long claimed /^-discrepancy for clusters (where -B spe c > P&t was claimed) 



has been resolved (Bahcall and Lubin 1994). The gas and galaxies trace each other both in 
their spatial density distribution and in their energies, as expected for a quasi- hydrostatic 
equilibrium. 
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Gas Mass Fraction The ratio of the mass of gas in clusters to the total virial cluster 
mass (within ~ 1.5/i -1 Mpc) is observed to be in the range 

^~0.03-0.15/i- L5 , (38) 

with a median value of 

<^P(< 1.5/T 1 Mpc)) median ~ 0.07k- 1 - 5 (39) 

(Jones and Forman 1992; White et al. 1993; White and Fabian 1995; Lubin et al. 1996). The 
implications of this result, which shows a high fraction of baryons in clusters, is discussed in 
§§ 

The total gas mass in clusters, ~ 1O 13 -1O 14 /i -2 ' 5 M , is generally larger than the total 
mass of the luminous parts of the galaxies (especially for low values of h) . With so much gas 
mass, it is most likely that a large fraction of the intracluster gas is of cosmological origin 
(rather than all the cluster gas being stripped out of galaxies). Additional optical-X-ray 
correlations of clusters are summarized in § |3.7| . 



3.6. Metal Abundance in Intracluster Gas 

The iron abundance in the intracluster gas is observed to be ~ 0.3 solar, with only small 
variations (± ~ 0.1) from cluster to cluster (e.g., Jones and Forman 1992). A strong corre- 
lation between the total iron mass in clusters and the total luminosity of the E + SO cluster 
galaxies is observed (Jones and Forman 1992). The metal enrichment of the intracluster gas 
is likely caused by gas stripped out of the elliptical galaxy members. 

The iron abundance profile as a function of radius from the cluster center is generally flat, 
i.e., a constant abundance at all radii (except for some poor, low-mass clusters dominated 
by a single massive galaxy). 

No evolution is observed in the overall iron abundance of clusters to z ~ 0.4 (Mushotzky 
1996). 

Recent results using ASCA observations (Mushotzky 1996) of different element abun- 
dances in nearby clusters (O, Ne, Mg, Si, S, Ar, Ca) suggest a SNII origin for the metals 
(resulting from early massive stars) rather than the expected SNIa. These new results will 
be expanded in the near future as additional accurate X-ray data become available and will 
provide further clues regarding the origin of the metallicity of the intracluster gas. 



-19- 



3.7. X-Ray— Optical Correlations of Cluster Properties 

Some observed correlations between X-ray and optical properties are listed in Table |7| 
(Bahcall 1977a,b; Edge and Stewart 1991; David et al. 1993; Lubin and Bahcall 1993). 

3.8. The X-Ray Luminosity Function of Clusters 

The X-ray luminosity function of clusters (the number density of clusters with X-ray 
luminosity L x to L x + dL x ) is approximately (Edge et al. 1990) 

$ x (L x )cLL x ~ 2.7 x 10- 7 (L X /10 44 )- L65 exp(-L x /8.1 x !0 u )(dL x /10 u ) Mpc -3 (40) 

for h = 0.5, where L x is the 2-10-keV X-ray luminosity in units of ergs -1 (for h = 0.5). The 
luminosity function can also be approximated as a power law (Edge et al. 1990) 

<5> x (L x )dL x ~ 2.2 x 10- 7 (L x /10 44 )- 2 - 17 (dL x /10 44 ) Mpc -3 (h = 0.5). (41) 

The number of X-ray clusters with X-ray luminosity brighter than L x is approximately 

n c (> L x ) ~ 2 x 10~ 7 (L X /10 44 )- L17 Mpc~ 3 (h = 0.5). (42) 

The observed evolution of the X-ray cluster luminosity function suggests fewer high-luminosity 
clusters in the past [z <: 0.5) (Edge et al. 1990, Henry et al. 1992). Additional data is re- 
quired, however, to confirm and better assess the cluster evolution. 

3.9. Cooling Flows in Clusters 

Cooling flows are common at the dense cores of rich clusters; X-ray images and spectra of 
~ 50% of clusters suggest that the gas is cooling rapidly at their centers (Sarazin 1986; Fabian 
1992). Typical inferred cooling rates are ~ lOOA^ /yr. The gas cools within r 100/i _1 kpc 
of the cluster center (generally centered on the brightest galaxy). The cooling flows often 
show evidence for optical line emission, blue stars, and in some cases evidence for colder 
material in HI or CO emission, or X-ray absorption. 

3.10. The Sunyaev-Zeldovich Effect in Clusters 

The Sunyaev-Zeldovich (1972) effect is a perturbation to the spectrum of the cosmic 
microwave background radiation as it passes through the hot dense intracluster gas. It is 
caused by inverse Compton scattering of the radiation by the electrons in the cluster gas. 
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Table 7: Correlations Between X-Ray and Optical 
Properties. a 

Properties Correlation 

a r -T a r (km s" 1 ) ~ (332 ± 52) [kT (kev)] a6±ai 

T-N . 5 kT (keV)~ 0.37V ° 5 95±ai8 

L x -N . 5 L x (bol)~ 1.4 x lO 40 ^ 16 ^- 15 ^- 2 

L x -f sp L x (bo\)^ 0.6 x 10 43 /- 2 - 16±an />- 2 

f sp -T f sp ~ 1.2[*T (kev)]-°- 94±0 - 38 

T-L x A;T (keV)~ O^L^bol^VlO 40 ] - 29 ™ 04 

a ay is the galaxy line-of-sight velocity disper- 
sion in the cluster (km s -1 ). T is the temper- 
ature of the intracluster gas [kT (keV)]. A/0.5 
is the central galaxy density in the cluster 
[number of galaxies brighter than m 3 + 2 m , 
within r < 0.5h^ = of the clus- 

ter center (Bahcall 1977a; Edge and Stew- 
art 1991)]. L x (bol) is the bolometric X-ray 
luminosity of the cluster (ergs -1 ). / sp is 
the fraction of spiral galaxies in the cluster 
1.5k' 1 Mpc) (Bahcall 1977b; Edge and 
Stewart 1991). Typical uncertainties of the 
coefficients are ~ 50% (see references). 
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At the long-wavelength side of the background radiation spectrum, the hot gas lowers 
the brightness temperature seen through the cluster center by the fractional decrement 

— = -2r ^, 43 

where T = 2.73 K is the microwave radiation temperature, To is the Thomson scattering 
optical depth through the cluster (r = <Jt f n e dl, where ot is the Thomson scattering cross 
section and dl is the distance along the line of sight), T x is the intracluster gas temperature, 
and m e is the electron mass. 

For typical observed rich cluster parameters of L x ~ 10 M h~ 2 ergs -1 , R c ~ 0.2/i _1 Mpc, 
and kT x ~ 4 keV, the bremsstrahlung relation (L x oc n 2 e B? c T® b , § |3.1|) implies a central 
gas density of n e ~ 3 x lO" 3 /^ 1 / 2 electrons cm -3 , thus yielding r — 3 x lO" 3 /^ 1 / 2 [r = 
. 64n e ( cm" 3 ) R c (kpc) ] . Therefore 

^~-6x lO' 5 ^ 1 / 2 . (44) 

This temperature decrement remains constant over the cluster core diameter 

2H R C 0.5 



cz z 

and decreases at larger separations. 



arcmin (45) 



The effect has been detected in observations of some rich, X-ray luminous clusters (e.g., 
Coma, A665, A2163, A2218, CI 0016+16) (Birkinshaw et al. 1984; Jones et al. 1993; 
Wilbanks et al. 1994; Herbig et al. 1995). 



4. The Baryon Fraction in Clusters 

Rich clusters of galaxies provide the best laboratory for studying the baryon fraction 
(i.e., the ratio of the mass in baryons to the total mass of the system) on relatively large 
scales of ~ Mpc. The mass of baryons in clusters is composed of at least two components: 
the hot intracluster gas (§|3|) and the luminous parts of the galaxies. 

The baryon fraction in clusters is therefore 

> M gas + M stars ^ Q Q7h -i. s + Q5 ^ (46) 



where the first term on the right-hand side represents the gas mass ratio (§ |3.5| ) and the 
second term corresponds to the stellar (luminous) contribution. 
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The baryon density required by big-bang nucleosynthesis is (Walker et al. 1991) 

fi 6 (BBN) ~ 0.015/T 2 . (47) 

Comparison of the above relations indicates that if Q m = 1 then there are many more baryons 
observed in clusters than allowed by nucleosynthesis. In fact, combining the two relations 
yields an Q m value for the mass-density of 

n ' m ~ 0.07h-™ + 0.05 ~ 007^ + ^05 ~ °' 2 (48) 

for the observed range of h ~ 0.5—0.8. Therefore, the baryon density given by nucleosynthesis 
and the high baryon content observed in clusters (mainly in the hot intracluster gas) suggest 
that Q m ~ 0.2. This assumes, as expected on this large scale (and as seen in simulations), 
that the baryons are not segregated relative to the dark matter in clusters (White et al. 
1993). 

Figure compares the observed gas mass fraction in clusters (oc baryon fraction) with 
expectations from cosmological simulations of Q m = 1 and Q m = 0.45 cold-dark-matter 
(CDM) flat models (Lubin et al. 1996). The results show, as expected, that the Q m = 1 
model predicts a much lower gas mass fraction than observed, by a factor of ~ 3. A low- 
density CDM model with Q m ~ 0.2—0.3 (in mass) best matches the data, as expected from 
the general analysis discussed above (see White et al. 1993; White and Fabian 1995; Lubin 
et al. 1996). 

In summary, the high baryon fraction observed in clusters suggests, independent of any 
specific model, that the mass density of the universe is low, Q m ~ 0.2—0.3. This provides a 
powerful constraint on high-density (fl m = 1) models; if Q m = 1, a resolution of this baryon 
problem needs to be found. 



5. Cluster Masses 

The masses of clusters of galaxies within a given radius, M(< r), can be determined 
from three independent methods: 

a) Optical: galaxy velocity dispersion (and distribution) assuming hydrostatic equilibrium 

(S); 

b) X-Rays: hot gas temperature (and distribution) assuming hydrostatic equilibrium 
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Fig. 2. — Observed and simulated gas mass fraction (M gas /M) vs. line-of-sight velocity 
dispersion for rich clusters (Lubin et al. 1996). A typical la uncertainty is shown. The 
simulated results (open squares) present the dark matter velocity dispersion; the galaxy 
velocity dispersion is lower by b v ~ 0.8 for SCDM and by b v ~ 0.9 for LCDM. (a) SCDM 
(fi m = l,h = 0.5); (b) LCDM (Q m = 0.45, h = 0.6). 
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c) Lensing: gravitational distortions of background galaxies (Tyson et al. 1990, 1996; Kaiser 
and Squires 1993). This method determines directly the surface mass overdensity. 

The first two methods were discussed in §§^ and |3|; the last (and newest) method is 
discussed in this book by Narayan. (See also Tyson et al. 1990, 1996; Kaiser and Squires 
1993; Smail et al. 1995, and references therein). 

The galaxy and hot gas methods yield cluster masses that are consistent with each 
other as discussed in §|3| (with o 2 ;, ~ kT/nm p ). Gravitational lensing results, which provide 
direct information about cluster masses, are available only for a small number of clusters 
so far (with data rapidly increasing). For all but one of the clusters the masses determined 
from lensing are consistent, within the uncertainties, with the masses determined from the 
galaxies and hot gas methods (see Bahcall 1995 for a summary). Some differences in masses 
between the different methods are expected for individual clusters due to anisotropic veloci- 
ties, cluster orientation (yielding larger lensing surface mass densities for clusters elongated 
in the line-of-sight, and vice versa), and sub-structure in clusters. On average, however, all 
three independent methods yield cluster masses that are consistent with each other. This 
triple check on cluster masses provides strong support for the important cluster mass deter- 
minations. 

The masses of rich clusters ranee from ~ 10 14 to ~ lO 15 /^ 1 M within 1.5/T 1 Mpc 
radius of the cluster center. When normalized by the cluster luminosity, a robust mass-to- 
light ratio is determined for nearby clusters, with only small variations from cluster to cluster 

M Mr, 

— (clusters) ~ 300 ± 100/i— - . (49) 
Lb L & 

This result is similar to the one obtained from the baryon fraction in §[|. 

If, as desired by theoretical arguments, the universe has critical mass density Q m = 1, 
than most of the mass in the universe cannotbe associated with galaxies, groups, and clusters; 
the mass distribution in this case would be strongly biased (i.e., mass does not follow light, 
with the mass distributed considerably more diffusely than the light). 



6. Where is the Dark Matter? 

A recent analysis of the mass-to-light ratio of galaxies, groups and clusters by Bahcall, 
Lubin and Dorman (1995) suggests that while the M/L ratio of galaxies increases with 
scale up to radii of R ~ 0.1-0.2/i _1 Mpc, due to the large dark halos around galaxies (see 
Fig. [3]; also Ostriker et al. 1974), this ratio appears to flatten and remain approximately 
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Fig. 3. — Mass-to-light ratio of spiral and elliptical galaxies as a function of scale (Bahcall, 
Lubin and Dorman 1995). The large boxes indicate the typical (~ la) range of M/L B for 
bright ellipticals and spirals at their luminous (Holmberg) radii. (Lb refers to total corrected 
blue luminosity; see text.) The best-fit M/L B oc R lines are shown. 
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constant for groups and rich clusters, to scales of ~ 1.5 Mpc, and possibly even to the larger 
scales of superclusters (Fig. f|). The flattening occurs at M/L B ~ 200-300/i, corresponding 
to Q m ~ 0.2. This observation suggests that most of the dark matter is associated with 
the dark halos of galaxies. Unlike previous expectations, this result implies that clusters 
do not contain a substantial amount of additional dark matter, other than that associated 
with (or torn-off from) the galaxy halos, and the hot intracluster medium. Bahcall et al. 
(1995) suggest that the relatively large M/L B ratio of clusters (~ 300/i) results mainly 
from a high M/Lb ration of elliptical/SO galaxies. They show (Fig. |3|) that ellipticals have 
an M/Lb ratio that is approximately 3 to 4 times larger than typical spirals at the same 
radius [{M/L B ) S ~ 100/i and (M/L B ) e ~ 400/i within r <, 200/i^ 1 Kpc]. Since clusters are 
dominated by elliptical and SO galaxies, a high M/Lb ratio results. 

Unless the distribution of matter is very different from the distribution of light, with 
large amounts of dark matter in the "voids" or on very large scales, the above results suggest 
that the mass density in the universe may be low, Q m ~ 0.2 (or Q m ~ 0.3 for a small bias 
of b ~ 1.5, where the bias factor b relates the overdensity in galaxies to the overdensity in 
mass: b = (Ap/ p) gal /(Ap/ p) m ). 

7. The Mass Function of Clusters 

The observed mass function (MF), n(> M), of clusters of galaxies, which describes the 
number density of clusters above a threshold mass M, can be used as a critical test of theories 
of structure formation in the universe. The richest, most massive clusters are thought to 
form from rare high peaks in the initial mass-density fluctuations; poorer clusters and groups 
form from smaller, more common fluctuations. Bahcall and Cen (1993) determined the MF 
of clusters of galaxies using both optical and X-ray observations of clusters. Their MF is 
presented in Figure |5|. The function is well fit by the analytic expression 

n(> M) = 4 x 10~ 5 (M/M*)- 1 exp(-M/M*)/i 3 Mpc" 3 , (50) 

with M* = (1.8 ± 0.3) x 10 14 /i _1 M , (where the mass M represents the cluster mass within 
1.5/T" 1 Mpc radius). 

The observed cluster mass function is compared in Figure |5] with expectations from 
different cold-dark-matter cosmologies using large-scale simulations (Bahcall and Cen 1992). 
The comparison shows that the cluster MF is indeed a powerful discriminant among models. 
The standard CDM model (Q m = 1) cannot reproduce the observed MF for any bias param- 
eter; when normalized to the COBE microwave background fluctuations on large scales, this 
model produces too many massive clusters, unseen by the observations. A low-density CDM 
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Fig. 4. — Composite mass-to-light ratio of different systems — galaxies, groups, clusters, and 
superclusters — as a function of scale (Bahcall et al. 1995). The best-fit M/Lb oc R lines 
for spirals and ellipticals (from Fig. are shown. We present median values at different 
scales for the large samples of galaxies, groups and clusters, as well as specific values for 
some individual galaxies, X-ray groups, and superclusters. Typical la uncertainties and la 
scatter around median values are shown. Also presented, for comparison, are the M/Lb (or 
equivalently Q) determinations from the cosmic virial theorem, the least action method, and 
the range of various reported results from the Virgocentric infall and large-scale bulk flows 
(assuming mass traces light). The M/Lb expected for Q = 1 and Q = 0.3 are indicated. 
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Fig. 5. — Cluster mass functions from observations and from CDM dimulations (Bahcall and 
Cen 1992). 
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model on the other hand, with Q m ~ 0.2—0.3 (with or without a cosmological constant), 
appears to fit well the observed cluster MF (Fig. ||). 

8. Quasar-Cluster Association 

Imaging and spectroscopic data (Yee et al. 1987; Ellingson et al. 1991; Yee et al. 1992) 
indicate that quasars are found in environments significantly richer than those of average 
galaxies. The data show a positive association of quasars with neighboring galaxies. 

Optically selected quasars to z 0.7 exhibit a quasar-galaxy cross-correlation function 
amplitude, A qg , that is approximately 2.3 times stronger than the galaxy-galaxy correlation 
amplitude (to separations r < 0.25/i _1 Mpc): 

(A*,) z> 2.3(A gg ) ~ 46. (51) 

This excess correlation suggests that the quasars are typically located in groups of galaxies 
with a mean richness 

rl.5 

(N R ) =n„ A qg r~ 1 - S A-Kr 2 dr ~ 12 galaxies (52) 
Jo 

(where n g ~ 0.015 Mpc -3 is the mean density of galaxies). The range of individual group 
richnesses is, however, wide. 

Radio-loud quasars at z < 0.5 are found in similar environments to those of the op- 
tical quasars above. At 0.5 < z < 0.7, the radio quasars appear to be located in richer 
environments, with 

(Agg) ~ 8(A gg ) ~ 160 (radio quasars, 0.5 <^ z <> 0.7). (53) 

This cross-correlation amplitude corresponds to a mean environment of rich clusters (R ~ 0, 
N R ~ 40). Radio quasars at these redshifts are thus typically found in rich clusters. 

The average galaxy velocity dispersion of the parent clusters associated with the quasars 
(Ellingson et al. 1991; Yee et al. 1992) is 

o r ~ 500 kms" 1 . (54) 

The observed auto-correlation function of optically selected quasars is approximately (Shaver 
1988) 

i qq {r, z ~ 0) ~ 10 2±a2 [r(Mpc)]- LS (55) 
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The quasar correlation strength is intermediate between the correlation of individual galaxies 
and the correlation of rich clusters. This correlation strength is consistent with the quasars lo- 
cation in groups of the above mean richness, as would be suggested by the richness- dependent 
cluster correlation function (§P~D|) . The quasars may thus trace the correlation function of 
their parent clusters (Bahcall and Chokshi 1991). 

Similar results are observed for the association of radio galaxies with groups and clusters. 
This association explains the observed increase in the strength of the radio galaxy correlation 
function over the general galaxy correlations (Bahcall and Chokshi 1992). 

9. Superclusters 

Rich clusters of galaxies are powerful tracers of the large-scale structure of the universe 
(Bahcall 1988, Peebles 1993). I summarize in the sections below the use of clusters in tracing 
the large-scale structure; I include superclusters (§[5|), statistics of the cluster correlation 
function (§PH), and peculiar motions on large scales (§|TT|) . 

9.1. Supercluster Properties 

Redshift surveys of galaxies reveal that superclusters are very large, high-density sys- 
tems of galaxies that are flattened or filamentary in shape, extending to tens of Mpc. The 
superclusters appear to surround large under-dense regions ("voids") of comparable sizes 
creating a "cellular-like" morphology of the universe on large scales (Gregory and Thomp- 
son 1978; Gregory et al. 1981; Chincarini et al. 1981; Giovanelli et al. 1986; de-Lapparent 
et al. 1986; da Costa et al. 1988; Rood 1988; Schectman et al. 1996; Landy et al. 1996). 

Large scale superclusters have been identified very effectively by rich clusters of galaxies 
(Bahcall and Soneira 1984), like high mountain peaks tracing mountain chains. Superclus- 
ters are generally defined as clusters of rich clusters of galaxies above a given spatial density 
enhancement /. Here / = n c (SC)/n c , where n c (SC) is the number density of clusters in 
a supercluster and n c is the mean number density of clusters. The observed superclusters 
are large flattened systems, extending to ~ 150/i -1 Mpc in size. The superclusters typi- 
cally contain several clusters. The high-density superclusters generally surround low-density 
regions. 

A complete catalog of superclusters was constructed by Bahcall and Soneira (1984) to 
z < 0.08. The catalog identifies all superclusters that have a spatial density enhancement 
/ > 20. The mean number density of the Bahcall-Soneira superclusters is ~ 10~ 6 h 3 Mpc -3 , 
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Table 8: Global Properties of Bahcall-Soneira Superclus- 
ters. 



Property 



/ = 20 superclusters 



Number density of SCs 
Number of clusters per SC 
Fraction of clusters in SCs 
Size of largest SC 
SC shape 

Volume of space occupied by SCs 



~ 10- 6 /i 3 Mpc- 3 
2-15 clusters 
54% 

~ 150/i" 1 Mpc 
Flattened 

-3% 



with a mean separation between superclusters of ~ 100/i 1 Mpc. A summary of the main 
properties of the superclusters is presented in Table |[ 

The superclusters trace well the structure observed in the more detailed, but smaller, 
galaxy redshift surveys. 



9.2. Superclusters and Pencil-Beam Surveys 

Observations of the redshift distribution of galaxies in narrow (~ 40 arcmin) pencil-beam 
surveys to z <, 0.3 (Broadhurst et al. 1990; hereafter BEKS) reveal a highly clumped and 
apparently periodic distribution of galaxies. The distribution features peaks of galaxy counts 
with an apparently regular separation of 128 Mpc, with few galaxies between the peaks. 
What is the origin of this clumpy, periodic distribution of galaxies? What does it imply for 
the nature of the large-scale structure and the properties discussed above? Bahcall (1991) 
investigated these questions observationally, by comparing the specific galaxy distribution 
with the distribution of known superclusters. 

Bahcall showed that the observed galaxy clumps originate from the tails of large super- 
clusters (§ |9.1|) . When the narrow-beams intersect these superclusters, which have a mean 
separation of ~ 100 Mpc, the BEKS galaxy distribution is reproduced. 

The redshift distribution of the superclusters is essentially identical to the galaxy redshift 
distribution, i.e., it reproduces the observed peaks in the BEKS survey, for z & 0.1. This 
indicates that the galaxy clumps observed in the pencil-beam survey originate from these 
superclusters as the beam crosses the superclusters' surface. The main superclusters that 
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contribute to the clumps were identified. For example, the first northern clump originates 
from the Coma-Hercules supercluster (= the Great- Wall); the second northern clump is 
mostly due to the large Corona Borealis supercluster. 

The narrow-beam survey of BEKS is directed toward the north and south galactic poles. 
Some of the Bahcall-Soneira superclusters coincident with the BEKS peaks are located at 
projected distances of up to ~ 50-100 Mpc from the poles. This suggests that the high- 
density supercluster regions are embedded in still larger halo surfaces, ~ 100 Mpc in size, 
and that these large structures surround large underdense regions. The observed number of 
clumps and their mean separation are consistent with the number density of superclusters 
and their average extent (§|9.1|). 

The narrow widths of the BEKS peaks are consistent with, and imply, flat superclus- 
ters. From simulations of superclusters and pencil-beams, Bahcall, Miller, and Udomprasert 
(1996) find that the observed peak- widths distribution is consistent with that expected of 
randomly placed superclusters with <, 15 Mpc width and ~ 150 Mpc extent. 

The Bahcall-Soneira superclusters may exhibit weak positive correlations on scales 
~ 100-150 Mpc (Bahcall and Burgett 1986). This implies that the superclusters, and thus 
their related galaxy clumps, are not randomly distributed but are located in some weakly 
correlated network of superclusters and voids, with typical mean separation of ~ 100 Mpc. 
This picture is consistent with statistical analyses of the BEKS distribution as well as with 
the observational data of large-scale structure. The apparent periodicity in the galaxy distri- 
bution suggested by BEKS is expected to be greatly reduced when pencil-beams in various 
directions are combined; the scale reflects the typical mean separation between large super- 
clusters, ~ 100-150/i -1 Mpc, but with large variations at different locations. 

10. The Cluster Correlation Function 

The correlation function of clusters of galaxies efficiently quantifies the large-scale struc- 
ture of the universe. Clusters are correlated in space more strongly than are individual 
galaxies, by an order of magnitude, and their correlation extends to considerably larger 
scales (~ 50/i -1 Mpc). The cluster correlation strength increases with richness (oc luminos- 
ity or mass) of the system from single galaxies to the richest clusters (Bahcall and Soneira 
1983; Bahcall and West 1992). The correlation strength also increases with the mean spatial 
separation of the clusters (Szalay and Schramm 1985; Bahcall and Burgett 1986; Bahcall and 
West 1992). This dependence results in a "universal" dimensionless cluster correlation func- 
tion; the cluster dimensionless correlation scale is constant for all clusters when normalized 
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by the mean cluster separation. 

Empirically, two general relations have been found (Bahcall and West 1992) for the 
correlation function of clusters of galaxies, £j = Ajr -18 : 

Ai<xN i} (56) 

A ~ (0.4d 4 ) L8 , (57) 

where Ai is the amplitude of the cluster correlation function, N{ is the richness of the galaxy 
clusters of type i (§[2.2[), and di is the mean separation of the clusters. Here di = n$ 3 , where 
Hi is the mean spatial number density of clusters of richness iVj (§ |2.3[ ) in a volume-limited, 
richness-limited complete sample. The first relation, Eq. flSq), states that the amplitude of 
the cluster correlation function increases with cluster richness, i.e., rich clusters are more 



strongly correlated than poorer clusters. The second relation, Eq. (p7|) , states that the 
amplitude of the cluster correlation function depends on the mean separation of clusters (or, 
equivalently, on their number density); the rarer, large mean separation richer clusters are 



more strongly correlated than the more numerous poorer clusters. Eqs. fl56|) and (|57j ) relate 
to each other through the richness function of clusters, i.e., the number density of clusters as a 
function of their richness. Equation (|57j ) describes a universal scale-invariant (dimensionless) 
correlation function with a correlation scale r oi = A\ 8 ~ OAdi (for 30 < di < 90/i" 1 Mpc). 

There are some conflicting statements in the literature about the precise values of the 
correlation amplitude, A^. Nearly all these contradictions are caused by not taking account 
of Eq. (|5"6D . When apples are compared to oranges, or the clustering of rich clusters is 
compared to the clustering of poorer clusters, differences are expected and observed. 

Figure |] clarifies the observational situation. The A^di) relation for groups and clusters 
of various richnesses is presented in the figure. The recent automated cluster surveys of APM 
(Dalton et al. 1992) and EDCC (Nichol et al. 1992) are consistent with the predictions 
of Eqs. (|56D and (|57D, as is the correlation function of X-ray selected ROSAT clusters of 
galaxies (Romer et al. 1994). Bahcall and Cen (1994) show that a flux-limited sample of 
X-ray selected clusters will exhibit a correlation scale that is smaller than that of a volume- 
limited, richness-limited sample of comparable apparent spatial density since the flux-limited 
sample contains poor groups nearby and only the richest clusters farther away. Using the 
richness-dependent cluster correlations of Eqs. ( |5T)| ) and Q5"T|), Bahcall and Cen (1994) find 
excellent agreement with the observed flux-limited X-ray cluster correlations of Romer et al. 
(1994). 

The strong correlation amplitude of galaxy clusters, and the large scales to which clusters 
are observed to be positively correlated (~ 50— 100/i -1 Mpc), complement and quantify the 
superclustering of galaxy clusters discussed in §0. Clusters of galaxies are strongly clustered 
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Fig. 6. — The univ ersal dimensionless cluster correlations: the dependence of correlation am- 
plitude on mean separation (Bahcall and West 1992). Data points include different samples 
and catalogs of clusters and groups, as well as X-ray-selected and cD clusters. Quasars and 
radio galaxies, as represented by their parent groups, are also included. The APM results 
are presented; they are consistent with the expected relation. 
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in superclusters of large scales (§^), consistent with the strong cluster correlations to these 
scales (pH|). 

This fundamental observed property of clusters of galaxies — the cluster correlation 
function — can be used to place strong constraints on cosmological models and the density 
parameter Q m by comparison with model expectations. Bahcall and Cen (1992) contrasted 
these cluster observations with standard and nonstandard CDM models using large N-body 
simulations (400/i _1 box, 10 7 ' 2 particles). They find that none of the standard Q m = 1 CDM 
models can fit consistently the strong cluster correlations. A low-density (O m ~ 0.2—0.3) 
CDM- type model (with or without a cosmological constant), however, provides a good fit 
to the cluster correlations (see Figs. 0-|9|) as well as to the observed cluster mass-function 
(§0, Fig. |5|). This is the first CDM- type model that is consistent with the high amplitude 
and large extent of the correlation function of the Abell, APM, and EDCC clusters. Such 
low-density models are also consistent with other observables as discussed in this paper. The 
Q m constraints of these cluster results are model dependent; a mixed hot + cold dark matter 
model, for example, with Q m = 1, is also consistent with these cluster data (see Primack's 
chapter in this book). 

The CDM results for clusters corresponding to the rich Abell clusters (richness class R > 
1) with d = 55/i _1 Mpc are presented in Figure |7| together with the observed correlations 
(Bahcall and Soneira 1983; Peacock and West 1992). The results indicate that the standard 
Q m = 1 CDM models are inconsistent with the observations; they cannot provide either the 
strong amplitude or the large scales (<: 50/i _1 Mpc) to which the cluster correlations are 
observed. Similar results are found for the APM and EDCC clusters. 

The low-density, low-bias model is consistent with the data; it reproduces both the 
strong amplitude and the large scale to which the cluster correlations are detected. 

The dependence of the observed cluster correlation on d was also tested in the simu- 
lations. The results are shown in Figure |8] for the low-density model. The dependence of 
correlation amplitude on mean separation is clearly seen in the simulations. To compare 
this result directly with observations, I plot in Figure [5] the dependence of the correlation 
scale, r D , on d for both the simulations and the observations. The low-density model agrees 
well with the observations, yielding r a « 0.4c?, as observed. The Q m = 1 model, while also 
showing an increase of r with d, yields considerably smaller correlation scales and a much 
slower increase of r (d). 

What causes the observed dependence on cluster richness [Eqs. (p6|- p^1 ? The depen- 
dence, seen both in the observations and in the simulations, is most likely caused by the 
statistics of rare peak events, which Kaiser (1984) suggested as an explanation of the ob- 
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0=0.25, h = 0.75, b=1.3 



R (h-'Mpc) 



Fig. 7. — Two-point correlation function of Abell R > 1 clusters, with mean separation 
55/i -1 Mpc, from observations and the CDM simulations (Bahcall and Cen 1992). 
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Fig. 8. — Model dependence of the cluster correlation function on mean separation d (CDM 
simulation: = 0.2, h = 0.5, b = 1) (from Bahcall and Cen 1992). 
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Fig. 9. — Correlation length as a function of cluster separation, from both observations and 
simulations (Bahcall and Cen 1992). 
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served strong increase of correlation amplitude from galaxies to rich clusters. The correlation 
function of rare peaks in a Gaussian field increases with their selection threshold. Since more 
massive clusters correspond to a higher threshold, implying rarer events and thus larger mean 
separation, Eq. ( p7|) results. A fractal distribution of galaxies and clusters would also produce 
Eq. (|57D (e.g., Szalay and Schramm 1985). 



11. Peculiar Motions of Clusters 

How is the mass distributed in the universe? Does it follow, on the average, the light 
distribution? To address this important question, peculiar motions on large scales are studied 
in order to directly trace the mass distribution. It is believed that the peculiar motions 
(motions relative to a pure Hubble expansion) are caused by the growth of cosmic structures 
due to gravity. A comparison of the mass-density distribution, as reconstructed from peculiar 
velocity data, with the light distribution (i.e., galaxies) provides information on how well the 
mass traces light (see chapter by Dekel, and 1994). The basic underlying relation between 
peculiar velocity and density is given by 

V-v = -Q° m 6 5 m = -Q° m 6 5 9 /b (58) 

where 5 m = (Ap/p) m is the mass overdensity, 5 g is the galaxy overdensity, and b = S g /S m 
is the bias parameter discussed in A formal analysis yields a measure of the parameter 
(3 = Q^/b. Other methods that place constraints on (3 include the anisotropy in the galaxy 
distribution in the redshift direction due to peculiar motions (see Strauss and Willick 1995 
for a review). 

Measuring peculiar motions is difficult. The motions are usually inferred with the aid 
of measured distances to galaxies or clusters that are obtained using some (moderately- 
reliable) distance- indicators (such as the Tully-Fisher or D n — a relations), and the measured 
galaxy redshift. The peculiar velocity v p is then determined from the difference between the 
measured redshift velocity, cz, and the measured Hubble velocity, vh, of the system (the 
latter obtained from the distance- indicator): v p = cz — vh- 

A summary of all measurements of (3 made so far is presented in Strauss and Willick 
(1995). The dispersion in the current measurements of (3 is very large; the various determina- 
tions range from f3 ~ 0.4 to ~ 1, implying, for b ~ 1, Q m ~ 0.2 to ~ 1. No strong conclusion 
can therefore be reached at present regarding the values of (3 or Q m . The larger and more 
accurate surveys currently underway, including high precision velocity measurements, will 
likely lead to the determination of (3 and possibly its decomposition into fl m and b (e.g., Cole 
et al. 1994). 
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Fig. 10. — Differential three-dimensional peculiar velocity distribution of rich clusters of 
galaxies for four cosmological models (Bahcall, Gramann and Cen 1994). 
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Clusters of galaxies can also serve as efficient tracers of the large-scale peculiar velocity 
field in the universe (Bahcall, Gramann and Cen 1994). Measurements of cluster peculiar 
velocities are likely to be more accurate than measurements of individual galaxies, since 
cluster distances can be determined by averaging a large number of cluster members as well 
as by using different distance indicators. Using large-scale cosmological simulations, Bahcall 
et al. (1994) find that clusters move reasonably fast in all the cosmological models stud- 
ied, tracing well the underlying matter velocity field on large scales. The clusters exhibit 
a Maxwellian distribution of peculiar velocities as expected from Gaussian initial density 
fluctuations. The model cluster 3-D velocity distribution, presented in Figure [10], typically 
peaks at v ~ 600 km s -1 and extends to high cluster velocities of ~ 2000 km s . The 
low-density CDM model exhibits lower velocities (Fig. [10]). Approximately 10% of all model 
rich clusters (1% for low-density CDM) move with v 10 3 km s _1 . A comparison of model 
expectation with recent, well calibrated cluster velocity data (Giovanelli et al. 1996) is 
presented in Figure [TT] (Bahcall and Oh 1996). The comparison between models and obser- 
vations suggests that the cluster velocity data is consistent with a low-density CDM model, 
and is inconsistent with a standard Q m = 1 CDM model, since no high velocity clusters are 
observed. 

Cen, Bahcall and Gramann (1994) determined the expected velocity correlation function 
of clusters in different cosmologies. They find that close cluster pairs, with separations 
r <^ 10/i 1 Mpc, exhibit strong attractive motions; the pairwise velocities depend sensitively 
on the model. The mean pairwise attractive cluster velocities on 5/i _1 Mpc scale ranges 
from ~ 1700 km s" 1 for Vt, m = 1 CDM to ~ 700 km s" 1 for Q m = 0.3 CDM. The cluster 



velocity correlation function, presented in Figure [12|, is negative on small scales — indicating 
large attractive velocities, and is positive on large scales, to ~ 200/i~ 1 Mpc — indicating 
significant bulk motions in the models. None of the models reproduce the very large bulk 
flow of clusters on 150/i _1 Mpc scale, v ~ 689 ± 178 km s -1 , recently reported by Lauer and 
Postman (1994). The bulk flow expected on this large scale is generally < 200 km s _1 for 
all the models studied (fi m = 1 and Q m ~ 0.3 CDM, and PBI). 



12. Some Unsolved Problems 

Considerable progress has been made over the last two decades in the study of clusters 
and superclusters of galaxies, as described in these lectures. However, many problems still 
remain open. I highlight some of the unsolved problems in this field that are likely to be 
solved in the coming decade. Currently planned large redshift surveys of galaxies and clusters 
such as the Sloan Digital Sky Survey and the 2dF survey, deep optical and X-ray surveys 
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Fig. 11. — Observed vs. model cluster peculiar velocity functions (from Bahcall and Oh 
1996). The Giovanelli and Haynes (1996) data are compared with model expectations con- 
volved with the observational errors. Note the absence of a high velocity tail in the observed 
cluster velocity function. 
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Fig. 12. — Velocity correlation function of rich (R > 1) clusters of galaxies for three models. 
Error bars indicate the la statistical uncertainties (from Cen et al. 1994). 
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using HST, Keck, ROSAT, ASCA, and AXAF, among other, should allow a considerable 
increase in our understanding of the nature and evolution of these fundamental systems. At 
the same time, state of the art cosmological simulations to be available in the next decade 
(e.g., Ostriker, this book) should greatly enhance our ability to compare the observations 
with detailed expectations from various cosmologies and hopefully narrow down the correct 
cosmological model of our universe. 

Here is a partial list of some of the interesting unsolved problems in the field of clusters 
and superclusters of galaxies. 

Clusters of Galaxies 

• What is the mass distribution and its extent in clusters of galaxies? Using gravitational 
lensing distortions, one can determine the mass density profile, p m (r), and total cluster 
mass, M(r), of clusters and compare it with the distribution of galaxies and gas for a 
large sample of clusters. 

• Does mass follow light on these scales? If not, what is the bias factor as a function of 
scale, b(r)7 

• What is the implied density parameter from clusters, Q m (r)? 

• What is the accurate baryon-fraction in clusters and groups of galaxies, as a function 
of scale, f2fe/f2(r)? 

• What is the origin of the hot intracluster gas and its metallicity? 

• What is the evolution of clusters in the optical and in X-rays? 

• What are the cosmological implications from studies of clusters? 



Superclusters 



• What is, quantitatively, the morphology of superclusters and large-scale structure (su- 
perclusters, filaments, and void network)? 

• What is the dependence of the superclustering properties on galaxy luminosity, surface 
brightness, type (E, S), and system (galaxies versus clusters)? 

• What are the peculiar motions in superclusters and on large scales? 
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• What is the mass, and mass distribution, in superclusters and on large scales? Does 
mass follow light? 

• What is in the "voids" ? 

• What is Q m on large scales? 

• What is the baryon fraction in superclusters? 

• What is the time evolution of superclusters? 

• What are the constraints made by the observed superclusters and large-scale structure 
on cosmology and galaxy formation models? 

I expect that many of these questions will be addressed and possibly solved in the 
coming decade. 
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